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ABSTRACT 



Since the launch of INTEGRAL in late 2002 we have 
monitored the Galactic microquasar GRS 1915+105 with 
long exposures (^100 ks) pointings. All the observa- 
tions have been conducted simultaneously with other in- 
struments, in particular RXTEand the Ryle Telescope, 
and in some cases with others (Spitzer, Nangay, GMRT, 
Suzaku,...). We report here the results of 3 observations 
performed simultaneously with INTEGRAL, RXTE, the 
Ryle and Nangay radio telescopes. These observations 
show the so-called ly and A classes of variability during 
which a high level of correlated X-ray and radio vari- 
ability is observed. We study the connection between 
the accretion processes seen in the X-ray behaviour, and 
the ejections seen in radio. By observing ejection during 
class A, we generalise the fact that the discrete ejections 
in GRS 1915+105 occur after sequences of X-ray hard 
dips/soft spikes, and identify the most likely trigger of the 
ejection through a spectral approach to our INTEGRAL 
data. We show that each ejection is very probably the 
result of the ejection of a Comptonising medium respon- 
sible for the hard X-ray emission seen above 15 keV with 
INTEGRAL. We study the rapid variability of the source, 
an find that the Low Frequency Quasi-Periodic Oscilla- 
tions are present during the X-ray dips and show variable 
frequency. The ubiquity of this behaviour, and the fol- 
lowing ejection may suggest a link between the QPO and 
the mechanism responsible for the ejection. 

Key words: accretion, accretion disks — black hole 
physics — stars: individual (GRS 1915+105) — X-rays: 
stars — X-rays: observations (INTEGRAL, RXTE). 



1. INTRODUCTION 



GRS 1915+105 is probably the most spectacular high en- 
ergy source of our Galaxy. An extensive review on it 
can be found in Fender & Belloni (2004). To summa- 
rize, GRS 1915+105 is a microquasar hosting a black 
hole (BH) of 14.0 ± 4.4 Mq (Harlaftis & Greiner 2004), 



it is one of the brightest X-ray sources in the sky and 
it is a source of superluminal ejection (Mirabel & Ro- 
driguez 1994), with true velocity of the jets > 0.9c. The 
source is also known to show a compact jet during its pe- 
riods of low steady levels of emission (e.g. Fuchs et al. 
2003). Multi-wavelength coverage from radio to X-ray 
has shown a clear but complex association between the 
soft X-rays and radio/IR behaviour Of particular rele- 
vance is the existence of radio QPO in the range 20^0 
min associated with the X-ray variations on the same time 
scale (e.g. Mirabel et al. 98). These so called "30-minute 
cycles" were interpreted as being due to small ejections 
of material from the system, and were found to correlate 
with the disc instability, as observed in the X-ray band. 
Extensive observations at X-ray energies with RXTE have 
allowed Belloni et al. (2000) to classify all the observa- 
tions into 12 separate classes (labeled with greek names), 
which could be interpreted as transitions between three 
basic states (A-B-C): a hard state and two softer states. 
These spectral changes are, in most classes, interpreted 
as reflecting the rapid disappearance of the inner portions 
of an accretion disc, followed by a slower refilling of the 
emptied region (Belloni et al. 1997). 
In the X-ray timing domain, GRS 1915+105 also shows 
interesting features, such as the presence of Low or 
High Frequency Quasi-Periodic Oscillations (LFQPO, 
HFQPO) whose presence is, as observed in other micro- 
quasars, tightly link the X-ray behaviour. LFQPOs with 
variable frequency during classes showing cycles have 
been reported. Correlations between the frequency and 
some of the spectral parameters have been pointed out 
(e.g. Markwardt et al. 1999, Rodriguez et al. 2002a,b). 
GRS 1915+105 is also one of the first BH systems in 
which the presence of HFQPOs has been observed, first 
at ~ 65 — 69 Hz, and up to ^ 170 Hz (Morgan, Remil- 
lard & Greiner, 1997; Belloni et al. 2006). Although in 
all cases the frequencies of all types of QPOs seem to 
indicate a link with the accretion disc, their presence is 
always associated to the presence of a hard X-ray tail in 
the energy spectra. 

The link between the accretion and ejection processes 
is, however, far from being understood and different 
kinds of model are proposed to explain all observation- 
nal facts also including the X-ray low (0.1-10 Hz) fre- 



Figure 1. Multiwavelength light curves of GRS 
1915+105 on October 17-18 2004 showing two inter- 
val of class V variability, a) Kyle (line) at 15 GHz, and 
Nangay (point) at 2. 7 GHz, b) JEM-X 3-13keV,c) ISGRI 
18-100 keV, d) RXTE 2-60 keV. 
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Figure 2. Multiwavelength light curves of GRS 
1915+105 on November 2004 during class A variability. 
The panels are the same as in Fig. Q] with the exception 
that no Nangay data are available. 
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quency QPOs. It should be added that until the launch 
of INTEGRAL all studies had been made below 20 keV 
with the PCA onboard RXTE, bringing, thus, only few 
constraints to the behaviour of the hard X-ray emitter 
(hereafter called corona). Since the launch of INTE- 
GRAL in late 2002 we have monitored the Galactic mi- 
croquasar GRS 1915+105 with long exposure (~ 100 ks) 
pointings. All the observations have been conducted si- 
multaneously with other instruments, in particular RXTE 
and the Ryle Telescope, and in some cases with others 
(Spitzer, Nan9ay, GMRT, Suzaku,...), with the aim to try 
to understand the physics of the accretion-ejection phe- 
nomena, including, for the first time, the behaviour of the 
source seen above 20 keV up to few hundred keV. Obser- 
vations of GRS 1915+105 performed during the first year 
of the campaign by INTEGRAL and RXTE have been re- 
ported elsewhere (e.g. Fuchs et al. 2003, Hannikainen et 
al. 2005, Rodriguez et al. 2004, Hannikainen et al. 2006 
in prep), while the INTEGRAL results on some of the 
sources in this field have been subject to several publi- 
cations (Hannikainen et al. 2004, Rodriguez et al. 2005, 
Rodriguez Shaw & Corbel 2006, Fritz et al. 2006). We 
report here the results obtained during observations per- 
formed during A02 and A03 showing sequences of X- 
ray hard dips/soft spikes (the cycles), followed by radio 
flares. The details of these analysis and extended discus- 
sions can be found in Rodriguez et al. (2006b). 



2. DESCRIPTION OF THE 3 OBSERVATIONS 



We focus in this paper on 3 observations respectively 
taken on October 17-18 2004 (Obs. 1), November 15- 
16 2004 (Obs. 2), and May 13-14 2005 (Obs. 3). The 
multiwavelength light curves are shown in Fig. [T|to|3] 
The JEM-X light curves show in all cases the occur- 



rences of soft X-ray dips of different duration, followed 
by a short spike marking the return to a high degree of 
soft X-ray emission and variability (hereafter cycle). In 
all the following we focus on the moments of cycle ac- 
tivity and more specifically on the intervals during which 
we have simultaneous data at radio wavelengths, and/or 
RXTE data for the timing analysis. This is particularly 
relevant for Obs.l and 3 during w hich trans itions be- 
tween different classes are seen (Rod riguez et a l. 2006b). 
We, therefore, avoid including in our analysis classes for 
which the radio behaviour is not known. These ob serva- 
tions are classified following iBelloni et al. I (l2000l) as v, 
A, and /3 type cycle activity for Obs. 1, 2 and 3, respec- 
tively. In all three cases we observe at least one radio 
flare, which is indi cative of a small ejection of material 
(iMirabel et alj[T998.) . after a cycle. The high flux shown 
by Nangay at 2.7 GHz during the class i/ observation sug- 
gests that the radio flares follow systematically each cy- 
cle, as reported by Klein- Wolt et al. (2002). The presence 
of radio oscillations is also known in class /3, while no 
radio/X -ray connection had eve r been observed during a 
class A (iKlein-Wolt et alJl2002h . 

Despite of the similarities at soft X-ray energies (i.e. the 
presence of cycles) and in the radio domain (the flares), 
and the connection between the 2 domains, large differ- 
ences remain between the 3 observations at the time of 
the X-ray radio connection. While in class ly the 18-100 
keV light curve shows a rather high degree of hard X-ray 
emission, and occurrences of dips correlated with the soft 
X-ray dips, the level of > 18 keV emission in class A is 
much lower (by a factor of ^2), and shows no occurrence 
of dips (Fig|2]i. In the class (3 observation, although the 
class obviously changes along the observation, the level 
of hard X-ray seems to remain roughly constant. The de- 
gree of variability, however, evolves and increases. Hard 
X-ray dips correlated to the softer ones appear. As in 
class 1/ the amplitude of the variation between the pre dip 
level and the dip level is lower than that seen at softer 



Figure 3. Multiwavelength light curves of GRS 
1915+105 on May 2005 with the appearance of class (3 
variability halfway through the observation. The panels 
are the same as in Fig. |2] 
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Figure 4. A section on the X-ray cycle followed by the 
radio flare from class X. The 4 intervals from which the 
spectra were accumulated ( see text) are indicated. From 
top to bottom, the panels respectively represent the Ryle 
light curve, the 3-13 keV JEM-X light curve, the 18-100 
keV ISGRI light curve, andthe 18-100/3-13 keV hardness 
ratio. 
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X-rays. This behaviour had also been pointed out by 
iRodriguez et al.. (,2002a) during a class a observation. 



3. X-RAY SPECTRAL ANALYSIS 



3.1. Selection of Good Time Intervals 

In order to study the spectral evolution of GRS 1915+105 
through the cycle and the possible origin of the radio 
ejection, we divided the cycles into different intervals 
from which JEM-X and ISGRI spectra were accumu- 
lated. Each cycle was divided in, at least, 3 intervals 
based on the soft X-ray count rate and the 3-13/18-100 
keV hardness ratio. The intervals are defined as follow; 

• Class I': each cycle is composed of 4 distinct fea- 
tures. The soft X-ray dip, having a hard spectrum 
(interval va), a short precursor spike (interval j^s), 
a short dip (interval I'c) and the main spike (interval 
I'd), the last 3 having soft spectra^ Given the obser- 
vation of radio oscillations dKlein- Wolt et al]f20()2h . 
i.e. the occurrence of radio flares after each cycle, 
the different intervals from all cycles were further 
accumulated together before the fitting to increase 
the statistics. 

• Class A: inspection of the cycle shows slight differ- 
ences here, with, in particular, occurrences of a pair 
of cycles at some moments, while at other intervals 
the cycles are isolated. Since no X-ray/radio con- 
nection had been observed in this class before, we 
extracted spectra from the unique cycle that is fol- 
lowed by an ejection. As in class ly this cycle can be 
divided into 4 intervals, (intervals A^i, Xb, Ac, Xd 
respectively represent the dip, a precursor spike, a 
short dip, and the main recovery spike (see FiglDi. 



• The spectral behaviour of class f3 has been 
extensively studied at s oft X-ray energy (e.g 
Markwardt et al. 1999; R odriguez et a l. 2002b), de- 
tailed an alysis of the INTEGRAL data will be pre- 
sented in [Rodriguez et al.l (l2006bh . Here 3 intervals 
can be distinguished, the dips, a precursor spike fol- 
lowed by another dip of longer duration than in the 
previous classes, and the major spike. As in the pre- 
vious classes, the first dip is h ard, while all remai n- 
ing intervals have soft spectra (iBelloni et al.ll2000h . 

Note the A, B, C, D labels for the different intervals are 
no t related to th e spect ral states (A, B and C) identified 
bv IBelloni et"aLl (l2000h . In order to be able to compare 
the different spectra, we fitted them all with the same 
model consisting of a thermal component (ezdisk) and 
a Comptonised one (comptt). Fig. |5] shows, as an ex- 
ample, the sequence of spectra from class A. 



3.2. Results 

In all cases the evolution of the source through the cy- 
cles is more prononced in the soft X-rays (Fig. [T] to 
|3]l. At first sight we could think that the variations are 
then caused by changes in the accretion disc rather than 
in the Comptonised component. When fitting the differ- 
ent spectra with physical models, however, it seems to be 
the contrary. In class i/ evolution through i/^ to vd trans- 
lates into an apparent approach of the accretion disc, with 
an increasing temperature (disc undetectable in i^a, and 
reaching ^ 2.1 keV with R,„ ~6 Rq in vu)- Interest- 
ingly between and i^c the disc seems to get closer (12 
to 6 Rg), and hotter (1.21 to 1.5 keV), or at least remains 
constant within the errors. We calculated the 3-50 keV 



Figure 5. The 4 spectra and best models superimposed from class A. The individual components are also shown. The 
disc component is represented by the dashed line, while the Comptonised component is represented by the solid line. The 
evolution through the dips is quite obvious, with in particular a spectacular drop in the Comptonised component between 
B and C. 





unabsorbed fluxes of the two spectral components in the 4 
intervals. While the flux from the disk increases through 
the whole sequence, that of the Comptonised component 
is not that regular. It first increases from va to vb before 
decreasing by a factor of 2.5 in vc (reaching 1.0 x 10~® 
erg cm~^ s^^ ), and slowly recovers in vd (1-5 x 10~^ 
erg cm~^ s^^ ). 

The same behaviour is also observed in class A. While the 
transition between As, Ac and Ad is marked by a rela- 
tively constant value of the temperature of the accretion 
disc (around 1 .8-2 keV) and value of the inner radius (5-7 
Rg), a spectacular evolution of the Comptonised compo- 
nent, in particular bewteen A b and Ac is observed (this 
is exemplified by Fig. |5]). The evolution of the 3-50 keV 
unabsorbed fluxes, follows the same pattern as in class v. 
During the dip the 3-50 keV disc flux represents ^ 30% 
of the 3-50 keV flux, its contribution increases greatly 
during the other intervals (it is always greater than 70% 
of the 3-50 keV flux). Between As and Ac, the 3-50 keV 
flux of the Comptonised component is reduced by a factor 



of - 11 (from 2.2 X 10"*^ erg cm"^ g-i to 0. 19 x 10"*^ erg 
cm"'^ s^^ ), when we leave all spectral parameters free to 
vary, and with a lower limit of 2.7 (if the disc temperature 
is frozen to the same value as in As). Again after Ac the 
Comptonised component slightly recovers (1.32 x 10^^ 
erg cm^^ s^^ ). 

The spectral behaviour of the source during class (3 has 
been extensively studied in the past. It is interesting to 
note that, as for the previous two classes, the disc temper- 
ature sl2wl};jncreasesduring the dip, even after the spike 
(e.g. HwanketaDIMl), while the inner radius remains 
small, indicating the disc is close to the co mpact object. 
The transition, after the spike, to state A ^Belloni et alj 
12000 ) the very soft spectral state of GRS 1915H-105, sug- 
gest that the spectacular changes occurring at the spike 
are related to the Comptonised component. 



Figure 6. Dynamical power spectrum (top) and 
RXTE/PCA light curve (bottom) on a portion of a cycle 
from class v. This cycle is one of the two for which a 
radio ejection is observed. 
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4. TIMING ANALYSIS: LOW FREQUENCY 
QPOS 



Although most of the cycles seen with RKTE were 
not simultaneous with the cycles followed by radio 
ejections, we produced dynamical power spectra of the 
cycles with the view to study whether or not the soft 
X-ray dips were associated with LFQPO of variable 
frequency, as been seen i n class 0, and a, i/ an d 9 (e.g. 
[Rodriguez et al.' 2002a'b'; Vadawale et alj|2003|). Fig. 
shows the dynamical power spectrum of class v together 
with the associated light curve of the unique RXTE cycle 
simultaneous with the presence of radio ejection. The 
presence of a LFQPO (the thick white line) of varying 
frequency is quite obvious during the X-ray dip. The 
LFQPO remain until the precusor spik e, and then dis- 
appear s. As observed in other classes ([Rodriguez et al.l 
l2002allbl: IVadawale et al.l l2003h . the frequency slowly 
increases during the dip, following the slow increase of 
the source count rate. 

The dip preceding the radio ejection during class A was 
unfortunately not covered by RXTE. Assuming the same 
pattern repeats during all dips, we extracted a dynamical 
power spectrum from the unique sequence of dip covered 
by RXTE. This is shown in Fig. [T] together with the 
associated light curve. 



As in the other classes showing soft X-ray dips, a strong 
QPO with variable frequency appears during the dip (Fig. 
|7]i. Here, however, the frequency seems less tightly corre- 
lated to the X-ray count rate, since while during the dips 
the count rate remains rather constant, the frequency of 
the QPO is less stable. 



Figure 7. Dynamical power spectrum (top) and 
RXTE/PCA light curve (bottom) on the unique 2 cycles 
from class A covered by RXTE. 
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5. DISCUSSION AND CONCLUSION 



We have presented observations of GRS 1915H-105 made 
in simultaneity with INTEGRAL, RXTE, the Ryle tele- 
scope and in one occasion the Nangay telescope. These 
observations belong to classes v, A, and /3, which have 
the characteristic of showing sequences of soft X-ray dips 
followed by high level of X-ray emission with high vari- 
ability, sequences we referred to as cycles. While in some 
of these classes ejections (observed in the radio domain) 
had been observed to f ollow to recovery to a high level 
of X-ray emission (e.g. iMirabel et al.]|1998l) . the obser- 
vation of a radio flare during the class A observation is, 
to our knowledge, the first ever reported. This allows us 
to generalise the fact that there is always an ejection fol- 
lowing a cycle, provide d the X-ray dip is long enough 
(iKlein-Wolt et al.l l2002h . The duration of the dip dur- 
ing class A is 500s, while it is as low as 330s during 
class 0. This duration is in agreement with the results of 
iKlein-Wolt eTaP (l2002h who observed that for an ejec- 
tion to occur a > 100s hard state must occur. 
Detailed inspection of the cycles shows that the se- 
quences are not as simple as it seems at first sight, but that 
in all cases a precursor spike occurs prior to the recovery 
to the high X-ray level. Our spectral analysis of the differ- 
ent portions of the cycles shows that the same evolution 
roughly occurs during the cycles from different classes. 
During the X-ray dip the source is in its hard state (or so- 
cafled C state for GRS 1915H-105 Befloni et al. (2000)). 
The precursor spikes correspond to a simultaneous soft- 
ening of the source spectrum and an increase in the flux 
of the two main spectral components, the disc and the 
corona. The short dip immediately following the precur- 
sor corresponds to the disappearance of the corona in the 
two cases we analysed in detai l. The same h ad been ob- 
served in the case of class (IChatvlll998l) . Given the 
observation of ejected material soon after, we conclude 
that: 



1) The ejected material is the coronal material 

2) The true moment of the ejection is the precursor 
spike 

It has to be noted that IChatyl (Il998h : IRodriguez et alJ 
(I2002bh : l\^awale et alj (|2003|) had come to similar con- 
clusions during other classes. Another point compatible 
with this scenario is that if we estimate the lag between 
the radio and the X-ray, we obtain 0.28h and 0.26h during 
class V, 0.3 Ih in class A, and 0.29 and 0.34h during class 
/3, therefore the delay is comparable in all classes. Inter- 
estingly the same conclusion had been drawn in another 
microquasar XTE J155 0— 564. During its 2000 outburst 
IRodriguez et al.l (l2003l) have shown that the discrete ejec- 
tion following the maximum of the X-ray emission was 
compatible with the ejection of the coronal medium. 
Although the physical properties of the dips are not 
exactly the same, some other similarities raise inter- 
esting questions. The timing analysis of the classes 
showing cycles all show the presence of a transient 
LFQPO appearing during the dip, and with a variable 
frequency which may be somehow correlated to the X - 
rav flux dSwank et al.l [19981 IRodriguez etall l2002allbl: 
IVadawaleetal.ll2003h . The presence of LFQPO during 
hard states is quite common, and may be a signature of 
the physics occurring during the accretion and the ejec- 
tion of matter, since a steady compact jet is usually asso- 
ciated with this state. In the case of the cycles, the situ- 
ation may be quite different. The hard state is short and 
no compact jet is observed. Instead, a discrete ejection 
seems to take place at the end of each cycle. The pres- 
ence of LFQPO during the dip may, again, indicate that 
the QPO phenomenom and the accretion-ejection physics 
are linked; 

Recently Tagg er et al.l (|2004|) have proposed a "magnetic 
flood" scenario in trying to explain the occurrence of dips 
and ejections, based on a magn etic instability, the Ac - 
cretion Ejection Instabihty (AEI, Tag ger & Pellatlll999h . 
This instability has the effect of transporting angular mo- 
mentum and energy from the inner region of the disc, 
and emitting them perpendicularly and directly into the 
corona. A strong observational signature of this AEI is 
the presence of LFQPOs. In this framework the X-ray dip 
would correspond to the appearance of the AEI (equiv- 
alent to a transition to a hard state with appearance of 
LFQPO), and the ejection could be d ue to a reconnectio n 
event in the inner region of the disc (iTagger et alJl2004 . 
An effect of the latter would be to blow the corona and 
would end observationally as an ejection. Although none 
of our results brings proof of such a scenario, and that 
other models exist, this scenario is compatible with our 
results. In particular with the generalisation of the ejec- 
tion during all classes with dips, the relative similarities 
between the dips of all classes, both at radio and X-ray 
energies, and the constant presence of LFQPO during the 
dips, we feel that this model is a very promising one. 
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